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qI,^". a picture for thermodynamics of the glassy state was introduced recently by us (Phys. Rev. 

' Lett. 79 (1997) 1317; 80 (1998) 5580). It starts by assuming that one extra parameter, the effective 

\ temperature, is needed to describe the glassy state. This approach connects responses of macroscopic 

observables to a field change with their temporal fluctuations, and with the fluctuation-dissipation 
relation, in a generalized, non-equilibrium way. Similar universal relations do not hold between 
energy fluctuations and the specific heat. 

In the present paper the underlying arguments are discussed in greater length. The main part of 
the paper involves details of the exact dynamical solution of two simple models introduced recently: 
uncoupled harmonic oscillators subject to parallel Monte Carlo dynamics, and independent spherical 
spins in a random field with such dynamics. At low temperature the relaxation time of both models 
SIh \ diverges as an Arrhenius law, which causes glassy behavior in typical situations. In the glassy regime 

Ch . we are able to verify the above mentioned relations for the thermodynamics of the glassy state. 

^ ' In the course of the analysis it is argued that stretched exponential behavior is not a fundamental 

property of the glassy state, though it may be useful for fitting in a limited parameter regime. 
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1. INTRODUCTION 



■ O ' Thermodynamics is an old but very powerful subject. It applies to a wide variety of systems ranging from ideal gases 



to crystals and black holes. Important contributions to its development were made by Carnot, Clausius, Kelvin, and 
Boltzmann. Equilibrium thermodynamics, better called "thermostatics", is a well understood subject, and applied 
every day in many fields of science. The work of Gibbs showed its tremendous generality via its relation to statistical 
physics (i.e. partition sums). We shall explain, however, that precisely this success has been a stand in the way for 
the systems our interest: glasses. 

Non-equilibrium thermodynamics for systems close to equilibrium was worked out in the first half of this century. 
\^ ] Typical applications are systems with heat flows, electrical currents, and chemical reactions. The basic assumption 
is the presence of local thermodynamical equilibrium, and the basic task is to calculate the entropy production. 
Important contributions to this field were made by de Donder, Prigogine, de Groot and Mazur. 

Non-equilibrium thermodynamics for systems far from equilibrium has long been a field of confusion. A typical 
application is window glass. Such a system is far from equilibrium: a cubic micron of glass is neither a crystal nor an 
^ ordinary undercooled liquid. It is an undercooled liquid that, in the glass formation process, has fallen out of its own 
\ metastable equilibrium. The glassy state is inherently a non-equilibrium state: a substance that is a glass in daily life 
^ i: (time scale of years) would behave as a liquid on geological time scales. If each 500 years a picture would be taken of 
a window glass, then the movie composed of these pictures would look very much like a movie of a soap film. 

Until our recent works on this field, the general consensus reached after more than half a century of research was: 
Q ] Thermodynamics does not work for glasses, because there is no equilibrium. Even before going into any detail, it is 
O ■ clear that this conclusion itself is confusing, because thermodynamics should also apply outside equilibrium. Inspired 
^ I by the success of Gibbsian theory, the whole non-equilibrium part of thermodynamics had been forgotten! The correct 
• formulation should of course have been: Equilibrium thermodynamics does not work for glasses, because there is no 
r^J ^ equilibrium, surely a less surprising and non-embarrassing statement. (This history shows once more how regretful it 

. is that equilibrium thermodynamics did not get known under its most proper name, "thermostatics".) 
" " ' The negative conclusion about the applicability of thermodynamics was mainly based on the failure to understand 
the Ehrenfest relations and the Prigogine-Defay ratio. It should be kept in mind that, so far, the approaches leaned 
very much on equilibrium ideas. Well known examples are the 1958 Gibbs-DiMarzio H and the 1965 Adam-Gibbs |] 
papers, while a 1981 paper by DiMarzio has title "Equilibrium theory of glasses" and a subsection "An equilibrium 
theory of glasses is absolutely necessary" In our opinion such approaches are not applicable, due to the inherent 
non-equilibrium character of the glassy state. In the course of the present work we shall encounter lots of instances 
where such approaches indeed fail to describe the physics. Notice, however, that this immediately rules out the by 
far most discussed model glass, namely the Gibbs-DiMarzio theory as a viable model for a realistic glass. For 
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instance, it would predict the original Ehrenfest relations to be always satisfied, in violence with experiments to be 
discussed. 

In our view the current lack of a thermodynamic description is quite unsatisfactory, since so many decades in time 
are involved, ranging from the microscopic sub-picosecond regime to, for silicate-rich glasses, almost the age of the 
solar system, thus covering more than 25 decades. Naively we expect that each decade has its own dynamics, basically 
independent of the other ones. We shall find support for this point in the models that we shall investigate below. 

Near the glass transition a glass forming liquid exhibits smeared discontinuities in quantities such as the heat 
capacity, the expansivity and the compressibility. This look similar to continuous phase transitions of the classical 
type, i.e., with specific heat exponent a — 0, even though the analogy is not perfect, because the smaller specific 
heat value occurs below the glass transition. It was then investigated whether the jumps satisfy the two Ehrenfest 
relations (the analogs for second order transitions of the Clausius-Clapeyron relation of a first order transition). As 
reviewed recently [Q, it was found that the first Ehrenfest relation, involving the jump in the compressibility, is 
always violated, while the second one, involving the jump in the specific heat, is usually satisfied, but not always. 
It has become fashionable to combine these two relations by introducing the so-called Prigogine-Defay ratio 11. For 
equilibrium transitions this quantity should be equal to unity, and it was generally expected that it cannot take values 
below unity. In glasses typical values are said to lie in the range 2 < 11 < 5, even though a very careful experiment 
on polystyrene led to 11 « 1 ||] . 

Our recent studies have radically changed the view point. We have realized that the first Ehrenfest relation is 
automatically satisfied, the only subtlety being its proper interpretation. We have also put forward that the Maxwell 
relation and the second Ehrenfest relation are modified in the glassy state, due to lack of equilibrium |6|. 

We have investigated the possibility that, within a yet unknown class of systems, the glassy state is described by 
one extra state variable. This is basically the age of the system, or the cooling rate under which the glass has been 
formed. We realized that in thermodynamics the effective temperature Tg is a more useful extra parameter JTot . 

This approach has led to a picture for thermodynamic relations between values of macroscopic observables ^ |^ 
Later it was extended to their fluctuations | [lO| ]. The picture also incorporates the so-called fluctuation-dissipation 
relation (FDR) , put forward in works by Horner and Cugliandolo and Kurchan jlj] ; for a review, see . 

This relation has become a central point in research on off-equilibrium systems. Our more general approach shows 
that the effective temperature that occurs in thermodynamics and the one that occurs in the fluctuation-dissipation 
relation are almost identical. 

In the course of our work we gained insights from analytical results combined with educated guessing on the p-spin 
interaction spin glass. Some initial studies had the purpose to flnd the physical meaning of the non-equilibrium 
replica free energy in spin glass models p^ ]. It has turned out that replica theory provides the two-temperature 
off-equilibrium free energy that we shall discuss in much more general context . The basic drawback of the p-spin 
model is that dynamics is not solved in the activated regime. For a model of directed polymers in a correlated random 
potential the situation is a little better, but so far it also lacks a complete solution in the activated regime 
Another model is the backgammon model, for which the dynamics at zero field has been partly solved |]l8| |l^]. One 
could couple the system to a particle bath, and the chemical potential would play the role of an external field. So far 
this case remains to be worked out. 

More promising is a model of independent harmonic oscillators with parallel Monte Carlo dynamics, introduced 
recently by Bonilla, Padilla and Ritort |^ . For this model the Hamiltonian and thus the statics is trivial. Nevertheless, 
the exactly solvable dynamics exhibits interesting glassy aspects. Since it has a too simple behavior in a field, we 
have recently studied a related simple non-mean-field model with trivial statics and interesting dynamics, namely the 
parallel Monte Carlo dynamics of independent spherical spins in a quenched random field pO(| . We expect that both 
models lie in the same class as the lattice-gas models with kinetic constraints of Kurchan, Peliti, and Sellitto |2^ [p2[ . 
The latter model can, however, not be solved analytically. 

In this work we shall give details underlying the picture proposed in reference [|o|. In section |] we shall recall 
this picture. In section ^ we shall derive dynamical equations for averages, correlations and responses in a model 
of uncoupled harmonic oscillators subject to Monte Carlo dynamics, introduced in reference In section ^ we 

analyze these equations in the non-equilibrium low temperature regime. In section |^ we analyze the closely related 
model of uncoupled spherical spins, introduced in (lO|. We close with a discussion and summary. 
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2.. THERMODYNAMIC PICTURE FOR A SYSTEM DESCRIBED BY AN EFFECTIVE 

TEMPERATURE 



A state that slowly relaxes to equilibrium is characterized by t, the elapsed time, sometimes called "age" or "waiting 
time". For glassy systems this is of special relevance. For experiments on spin glasses it is known that non-trivial 
cooling or heating trajectories can be described by an effective age Yet we do not wish to discuss spin glasses 
in this work. They have an infinity of long time scales, or infinite order replica symmetry breaking. Their phase 
transition is continuous, and involves power laws. 

We shall restrict our treatment to systems with one diverging time scale, having, in the mean field limit, one step 
of replica symmetry breaking. They are systems with first-order- type phase transitions, with discontinuous order 
parameter, though usually there is no latent heat. However, the same approach applies to true first order glassy 
transitions that do have a latent heat. This occurs, for instance in the transition from low density amorphous ice to 
high density amorphous ice ||2^ |^^. Theoretically such behavior occurs in spin glasses in a transverse field, see e.g. 

@ [0 n n. 

We shall consider glassy transitions for liquids as well as for random magnets. The results map onto each other by 
interchanging volume V, pressure p, compressibility k — —dlnV/ dp, and expansivity a — dlnV/ dT, by magnetization 
A/, field H, susceptibility x = {l/N)dM/dH, and "magnetizabihty" a = {-l/N)dM/dT, respectively. 

The picture to be investigated in this work starts by describing a non-equilibrium state characterized by three 
parameters, namely T,H and the age effective temperature Te{t). As we shall see below, Te{t) enters naturally in the 
dynamical solution of the problem. For a set of smoothly related cooling experiments Ti{t) at different fields Hi, one 
may express the effective temperature as a continuous function: Tf,^i{t) Te{T, H). This sets a surface in (T, Tg, H) 
space, that becomes multi-valued if one first cools, and then heats. For covering the whole space one needs to do 
very many experiments, e.g., at different fields and at different cooling rates. The results should agree with findings 
from heating experiments and aging experiments. Thermodynamics amounts to giving differential relations between 
observables at nearby points in this space. 

For thermodynamics of glassy systems in the absence of currents, all previous results can be summarized by 
expressing the change in heat as 

dQ^TdS,p + T,dI (2.1) 

where Sep is the entropy of the fast or equilibrium processes (/3-processes) and I the configurational entropy of the 
slow or "configurational" processes (a-processes). This object is also known as information entropy or complexity. 
Both Sep and 2 are state functions in the sense that they depend on T, Tf., and on H or p. In particular, they are 
defined for any Tg, and, within the present framework of one effective parameter, they do not depend on the path 
along which this value was reached. 



Notice that our separation in eq. (2.1) goes according to time scales. In the common use of the word, the 
configurational entropy Sc is the entropy of the glass minus the entropy of the vibrational modes of the crystal . 
For polymers, in particular, it still includes short-distance rearrangements, which is a relatively fast mode. For the 
Gibbs-DiMarzio model it was confirmed numerically that Sc indeed does not vanish at any temperature, thus violating 
the Adam-Gibbs relation r ~ exp(const. /S'c) between time scale and configurational entropy Our I, on the other 
hand, only contains the slow components; the fast ones are supposed to be in equilibrium, and are counted in Sep- The 
properly formulated Adam-Gibbs relation should only refer to slow quantities, so it should read: r ~ exp(const/X). 
Its applicability remains an open issue. In a certain model glass with non-trivial fast and slow modes that has a 
Kauzmann transition it is actually satisfied |p| . 



In the presence of currents eq. (|2.l| ) would become dQ < TdS'cp -I- TedX. This decomposition is based on a system 
consisting of two parts, with a slow exchange of heat between them, so having two time scales. A well known case 
is a cup of coffee at temperature Te in a room at temperature T. In that case T is the entropy of the cup and the 
coffee, Sep the entropy of the air and matter in the room, and Q the heat of the combined system. To mention one 
case, cooling of the coffee in an isolated room will be described by dQ = and TgdJ = —TdSep < 0. 

It is both surprising and satisfactory that a glass can be described by the same general law. If also an effective 
pressure or field would be needed, then dQ is expected to keep the same form, but dW would change from its standard 
value —pdV for liquids, or —MdH for magnets. In the latter case it would become —MidH — M2dHe, where is 
the effective field, and Mi and M2 add up to M. Such a complication could be needed in a larger class of systems. 
It would make the picture technically a bit more difficult, and is the subject of current research. 
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1. First and second law 



For a glass forming liquid the first law dll = dQ + dW becomes 



dU = Td^cp + TedX - pdV 



(2.2) 



One can define the free enthalpy 



G = U- TS,p - T,I 



pV 



(2.3) 



that satisfies 



dG 



-SepdT - IdTe 



Vdp 



(2.4) 



The total entropy is 



(2.5) 



(We should stress that the total entropy is not equal to Sop + Tel/T; there are many reasons why this unsym- 
metric form is incorrect. Let us mention that if the probability distribution decomposes into fast and slow pro- 
cesses as P(fast, slow) = P(fast|slow)P(slow), then the standard expression S = — trPlnP leads to (2^) with 
Sep — trsiowP(slow)[— trtast-P(fast|slow) lnP(fast|slow)], just the entropy of the fast processes, averaged over the slow 
ones, and I = — trsiowP(slow) InP(slow), just the entropy of the slow processes.) 
The second law requires dQ < TdS, so 

(Te - T)dJ < (2.6) 
Since Tg = Te{T,p), and both entropies are functions of T, T^, and p, the expression (^.l[) yields a specific heat 
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dT 
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die 
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dn 



die 
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(2.7) 



In the glass transition region it holds that Tg « T. Since the derivatives of S'op and X are smooth functions, all factors, 
except dxTe, are basically constant. This leads to 



Cp — Ci + C; 



dTe 
dT 



(2., 



Precisely this form has been assumed half a century ago by Tool pCf] as starting point for the study of caloric behavior 
in the glass formati on r egion, and has often been used for the explanation of experiments [^T| js^. It is thus a direct 
consequence of eq. (2T). Let us mention that Tool uses the term "fictive temperature" for T,. 
For magnetic systems the first law brings 



dU = TdScp + TedJ - MdH 



(2.9) 



One can define the free energy 



F = U- TS,p - Tel 



(2.10) 



that satisfies 



dP 



'S'cpdP - XdTe 



MdH 



(2.11) 
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2. Modified Maxwell relation 



For a smooth sequence of cooling procedures of a glassy liquid, eq. ( |2.2| ) implies a modified Maxwell relation 
between macroscopic observables such as U{t,p) U{T,p) — U{T,Te{T,p),p) and V. This solely occurs since Te is 
a non-trivial function of T,p for the smooth set of experiments under consideration. 

The consistency relation d'^G/dTdp — d'^G/dpdT yields 

(2.12) 
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p dp 



Notice that difference relations as eq. (2.9), and the Legendre transformation that leads to (2.11), do not invoke the 
functional dependence Te{T,p), since they hold for any functional dependence, and even in absence of it. However, it 
does be com e relevant when dividing these equations by dT or dp, as was done to derive (2.12). 
Eq. (U) implies 



j.dS,p 



Eliminating dSep/dp leads to 

dU 
dp 



dV 
dp 



dp 



T dT 



dU 








dV 


T dp 
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^ dp 



T 



dl 
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p dp 
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dl 
dp 



T dT 



T 



dl 
dp 



This is the modified Maxwell relation between observables U and V. In equilibrium T, 
vanishes, and the standard form is recovered. 
Similarly, one finds for a glassy magnet 



dU 
dH 



M ~T 



dM 
'dT 



T, 



H 



dl 
'dH 



T 



[dH 



dl 
Tdf 



H 



dTe 

dT 



dl 
HdH 



(2.13) 
(2.14) 

T, so the right hand side 
(2.15) 



3. Modified Clausius-Clapeyron relation 

Let us consider a first order transition between two glassy phases A and B. An example could be the transition 
from low-density-amorphous ice to high-density-amorphous ice |2^. For the standard Clausius-Clapeyron relation 
one uses that the free enthalpy G is continuous along the first order phase transition line Pg{T). Since T, ^ T, it 
is actually not obvious that G should still be continuous there. The so far always confirmed fact that in mean field 
mod els replica theory brings the relevant physical free energy, leads us to expect that the generalized free enthalpy 



( p.3D is indeed continuous. 

Let us consider a first order transition between phases A and B, that have their own Tg, Sop and X. Let us denote 
the discontinuities in observables O of the two states as 

AO{T,pg{T)) = OA-OB (2.16) 
Taking O — G and differentiating AG = one gets 



dp 



j]]^^AS. • ^^^^^ 



T J J dT \dT 
5ep can be eliminated by means of (|2.3D . Using again that AG ~ 0, this yields 



I] (2.17) 

p 



, dpn AU+p„AV ^ fdT„ T„\ 

AV-^ = ^ hA — -I -l] (2.18) 

dT T \dT T J ^ ' 

where d/dT — d/dT + {dpcJdT)d/dp is the "total" derivative, i.e., the derivative along the transition line. This is 
the modified Clausius-Clapeyron relation. It would be very interesting to test this relation for ice. For that substance 
Mishima and Stanley | p5[ | have presented a thermodynamic construction of the free enthalpy or Gibbs potential G. 
It is, however, based on equilibrium ideas and does not involve the effective temperature in the amorphous phases. 
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In particular, it assumes the validity of the original Clausius-Clapeyron relation. We feel that the results are not 
the physically relevant ones, and that the analysis should be redone within our non-equilibrium thermodynamic 
framework. 

When phase A is an equilibrium undercooled liquid, and phase B is a glass, it holds that — T in phase A, and 
its X-terms will cancel from (2.18), so this relation reduces to 

dp^^AU + pAV T^ ,m 

dT T dT' ^ ' 

where Tg and T are properties of the glassy phase B. Notice that eq. (7) of ref. |6) contains a misprint in the prefactor 
of I. 

For standard glass forming liquids, there are no discontinuities in U and V. It then holds that along the glass 
transition line Te{T,pg{T)) — T, implying dT^/dT = 1, which indeed removes the I terms from the last two relations. 



4. Ehrenfest relations and Prigogine-Defay ratio 

In the glass transition region a glass forming liquid exhibits smeared jumps in the specific heat Cp, the expansivity 
a and the compressibility k. If one forgets about the smearing, one may consider them as true discontinuities, yielding 
an analogy with continuous phase transitions of the classical type. 

Following Ehrenfest one may take the derivative of /Sy{T,pg{T)) ~ 0. Using the definitions of a and k, given 
above, the result for a glass forming liquid may be written as 

Aa = Ak^ (2.20) 

dT ^ ' 



while for a glassy magnet 



^ dT ^ ' 



The conclusion drawn from half a century of research on glass forming liquids is that this relation is never satisfied 
pl| ]3^ ] Q|. This has very much hindered progress on a thermodynamical approach. However, from a theoretical 
viewpoint it is hard to imagine that something could go wrong when just taking a derivative. We have pointed out 
that this relation is indeed satisfied automatically but it is important say what is meant by k in the glassy state. 

Let us make an analogy with spin glasses. In mean field theory they have infinite order replica symmetry breaking. 
From the early measurements of Canella and Mydosh on AuFe it is known that the susceptibility depends loga- 
rithmically on the frequency, so on the time scale. The short-time value, called Zero-Field-Cooled (ZFC) susceptibility 
is a lower bound, while the long time value, called Field-Cooled (FC) susceptibility is an upper bound. Let us use the 
term "glassy magnets" for spin glasses with one step of replica symmetry breaking. They are relevant for comparison 
with glass forming liquids. For them the situation is worse, as the ZFC value is discontinuous immediately below 
Tg. { At H = one has xzfc = /3(1 — Qea), while xfc — /5(1 ^ (1 ~ '^\)(Iea) matches xpm = /3 at xi = 1.) This 
occurs since giving the system more time to react on the field, will lead to a really larger response, and it explains 
why already directly below the glass transition different measurements yield different values for k. These notions are 
displayed in figure 2.1 
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FIG. 2.1. Schematic plot of the field-cooled (FC) and zero-field-cooled (ZFC) susceptibility in realistic spin glasses and 
in glassy magnets, as function of temperature, in arbitrary units. In realistic spin glasses the infinite time or field-cooled 
susceptibility is larger than the short time or zero-field-cooled susceptibility. In magnetic analogs of realistic glasses the short 
time susceptibility even has a smeared discontinuity at the glass transition. In glass forming liquids the same happens for the 
compressibility. 

Previous claims about the violation of the first Ehrenfest relation can be traced back to the equilibrium thermo- 
dynamical idea that there is one, ideal k, to be inserted in ( ^.20| ). Indeed, investigators usually considered coohng 
curves V{T,pi) at a set of pressures pi to determine Aa and dpg/dT. (An alternative route, often followed in polymer 
physics, and leading to a very similar problem, is to change p at many constant values of T; then k depends strongly on 
the rate of change of p.) However, Ak was always determined in another way, often from measurements of the speed of 
sound, or by making more complicated pressure steps [|j . In equilibrium such alternative determinations would yield 
the same outcome. In glasses this is not the case: the speed of sound is a short-time process, and additional pressure 
steps modify the glassy state. Therefore alternative procedures are not allowed, and only the cooling curves V{T,pi) 
should be used. They constitute a liquid surface V\ic^nid(T , p) and a glass surface V^iass(2^iP) in {T,p,V) space. These 
surfaces intersect, and the first Ehrenfest relation is no more than a mathematical identity about the intersection line 
of these surfaces. It is therefore automatically satisfied The most careful data we came across were collected by 
Rehage and Oels for atactic polystyrene Q. In figure 2.2 we present those data in a 3-d plot, underlining our point 
of view. 

After submitting the original version of this paper, we realized that McKenna has stressed that in experiments on 
glasses the isothermal compressibility differs from the isochoral compressibility |^ . He also concludes that alternative 
experiments are not allowed, and that the first Ehrenfest relation indeed is merely a tautology. 
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FIG. 2.2. Data of the glass transition for cooling atactic polystyrene at rate 18 K/h, scanned from the paper of Rehage 
and Oels (1976): specific volume V {crn?/g) versus temperature T (K) at various pressures p (kbar). As confirmed by a 
polynomial fit, the data in the liquid essentially lie on a smooth surface, and so do the data in the glass. The first Ehrenfest 
relation describes no more than the intersection of these surfaces, and is therefore automatically satisfied. The values for the 
compressibility derived in this manner will generally differ from results obtained via other procedures. 

The second Ehrenfest relation derives from differentiating AU{T,pg{T)) = 0. The obtained relation will also be 
satisfied automatically. However, one then eliminates dU/dp by means of the Maxwell relation. In equilibrium this 
would yield 



ACp 

TgV 



Aa 



dpg_ 
dT 



(2.22) 



We have already discussed that outside equilibrium it is modified, see eq. ( ^.14 ). We thus obtain instead 
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die 

dT 
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dT 



Aa^ 
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die 
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df 
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(2.23) 



where dX/dT is the "total" derivative of the configurational entropy along the glass transition line. The last term is 
new and vanishes only at equilibrium. For magnets one gets 



AC 

iVT 



Aa— — f 
dT 



dT 



H 



dl_ 
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dHg dl 
H ^ 'dT 'dH 



(2.24) 



Along the glassy transition line the equality Te{T, Hg{T)) — T implies 
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die 

dH 



T dT 
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Combining the two original Ehrenfest relations one may eliminate the slope of the transition line, 
consider the so-called Prigogine-Defay ratio 



(2.25) 



This leads to 



n 



ACpAn 
TV{Aaf 



(2.26) 



For equilibrium transitions it should be equal to unity. Assuming that at the glass transition a number of unspecified 
parameters undergo a phase transition, Davies and Jones showed that 11 > 1 while DiMarzio showed that in that 
case the correct value is 11 = 1 In glasses typical experimental values are reported in the range 2 < 11 < 5. It 

was therefore generally expected that 11 > 1 is a strict inequality. 

We have pointed out, however, that, as the first Ehrenfest relation is satisfied but the second is not, it holds that 



n 



TVAa{dpg/dT) 



1 



1 



VAa 



1 



dT 



dX 
Ap 



(2.27) 



Depending on the set of experiments to be chosen, Apg/AT can be small or large, and 11 can also be below unity. 
Rehage and Oels found 11 — 1.09 k \ aX p — 1 kbar, using a short-time value for k Reanalyzing their data we 
find from ( ^.27 ), where the physically relevant k has been inserted, a value 11 = 0.77, which is, surprisingly enough, 
below unity ||6|]. 

The definition (2.26) of 11 looks like a combination of equilibrium quantities. This is misleading, however, since 



'igiass depends sensitively on how the experiment is done. We conclude that the commonly accepted inequality 11 > 1 
is based on equilibrium assumptions. Our theoretical arguments and the Rehage-Oels data show that such idea's are 
incorrect. In particular this rules out the Gibbs-DiMarzio model as a principally correct model for the glassy state. 
It is an equilibrium model, and as such it will e.g. lead to 11 = 1, in contradiction to experiments. 



5. Fluctuation formula 



The basic result of statistical physics is that it relates fluctuations in macroscopic variables to response of their 
averages to changes in external field or temperature. We have wondered whether such relations generalize to the glassy 
state. We have found arguments in favor of such a possibility both from the fluctuation-dissipation relation and by 
exactly solving the dynamics of model systems ]To[ |. Susceptibilities appear to have a non-trivial decomposition, that 
looks as being very general. Here we give arguments leading to it. 

Later we shall consider models where two fields Ha {a = 1,2) are present, and two magnetizations Ma (a = 1,2). 
In cooling experiments at fixed field H = {Hi,H2) it holds that Ma —Ma{T{t),Te{t, H), H). For thermodynamics 
one eliminates time to express Te{t,H) — > Te{T,H), implying Ma —Ma{T,Te{T,H),H). One may then expect three 
terms: 



Xab 
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N dHh 



fluct 
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it) 



it) 



(2.28) 



The first two are defined by 



xT'(t) + xTit)^ 



1 dMa 



N dHb 



T,T, 



(2.29) 



To calculate them separately, we switch from a cooling experiment to an aging experiment at the considered T, Tf. and 
-ff, by keeping, in Gedanken, T fixed from then on. The system will continue to age, expressed by Tg = Te{t]T, H). 
We may then use the equality 



dMa 




dMa 




dMa 


dT, 


dHb 


T,t 


dHb 




dT, 


T,H dHb 



T,t 



(2.30) 



We have conjectured |1C(] that the left hand side may be written as the sum of fluctuation terms for fast and slow 
processes, 



xlT\t) = 



1 dMa 
N dHb 



{6Ma{t)6Mb{t))f,,t , {6Ma{t)SMb{t)) 

slo 



T,t 



NT{t) 



NT,{t) 



(2.31) 
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The first term is just the standard equihbrium expression for the fast equihbrium processes. Notice that slow processes 
enter with their own temperature, the effective temperature. This decomposition is confirmed by use of the fluctuation- 
dissipation relation in the form to be discussed below. Combination of ( ^.29 ), (2.30) and (2.31) yields 



1 dMa 



N dT, 



TM dHb 



T,t 



(2.32) 



The fluctuation terms are instantaneous, and thus the same for aging and cooling. The loss term is a correction, 
related to an aging experiment. It measures the decrease of fluctuations below the glass transition, which will be 
small in the models to be discussed later on. 

In the models to be considered below, dynamics in the glassy phase is essentially independent of the actual T, 
leaving almost no difference between cooling and aging. This is due to the simplicity of the model. 

Since Tg ^ T, there occurs in eq. (2.28) also a new, conflgurational term 



1 dMa 



N dT, 



T.H dHh 



(2.33) 



It originates from the difference in the system's structure for cooling experiments at nearby flelds. For glass forming 
liquids such a term occurs in the compressibility. Its existence was anticipated in some earlier works. Goldstein |33| ] 
points out that V^iass depends stronger on the pressure of formation pform than on the one remaining after partial 
release of pressure, pfinai- Jackie |3^ then assumes that for infinitely slow cooling pform is the only additional system 
parameter, and argues that Akt — > Ak — /S.kt ^ dhiV / dpfo^^ — AadTg/dp and that this implies 11 = Akt/Ak > 1. 
He thus also considers one extra system variable, and also argues the existence of a conflgurational term. We do 
not wish to restrict to adiabatically slow cooling, and we do not agree with his conclusion on 11. Notice that our 
approach allows, in principle, to find the configurational term (2.32) for typical cooling procedures from construction 
of V{T,Te,p) in fuU (T, re,p)-space. 

From the analysis to be given below, we find no reason why such universal quasi-equilibrium relations could also 
hold between the specific heat and the energy fluctuations. In the models of the present paper the energy fluctuation 
arc smaller by one order of magnitude, and model dependent. The absence of such a general relation allowed us to 
apply the very same two-temperature approach to black holes, without obtaining a contradiction with their negative 
specific heat [BSl. 



6. Fluctuation-dissipation relation 

Nowadays quite some attention is payed to the fiuctuation-dissipation relation in the aging regime of glassy systems. 
It was first put forward in works by Horner [|l2| ||l^ and then by Cugliandolo and Kurchan iQ. This relation 
has become a central point in research on off-equilibrium systems, for a review, see ||l5| . 

Our formulation is that in the aging regime there holds the following relation between the cross correla- 
tion Cab{t,t') = {Oa{t)Ob{t')) — {Oa{t))(Ob{t')) of macroscopic observables Oait) and Ob{t'), and the response 
Gab{t,t') = S{Oa{t))/SHb{t') of Oa{t) to 'A short, small field change applied at an earlier time t': 

^^^^^Ml^f,^t')Gab{t,t') (2.34) 

with Te{t) being the effective temperature for the FDR, while in the equilibrium or short-time regime, T replaces 
Te- This relation has been confirmed numerically e.g. for a soft sphere glass [^. It is remarkable that the Te is a 
function of one of the times only. However, one should keep in mind that C and G typically have a t' /t scaling, while 
Te typically is a smooth function of Ini, a variable that basically equals Int' . 

One expects that Te{t) is close to the "thermodynamic" effective temperature Te{t). Let us show how this comes 
about. 



The two fiuctuation terms in eq. (2.31) are consistent with (2.34). To prove this, let us neglect switching effects 



(see section 3.4), and use the definition 

1 dMa 



N dHb 



= / Gab(t,t')dt' (2.35) 
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We split the integral up in the regions {t — Tp, t) and (0, t — Tp), where Tp is the time after which the fast or /? processes 
have died out. Their contribution has the equilibrium form, while in the second interval we may insert (2.34), which 
yields 



1 dM,, 



N dHb 



{SMjt)SMt{t)) 



fast 



NT{t) 



fe{t') dt' 



(2.36) 



We perform a partial integration, and can neglect the value at the lower boundary t — Q. In the remaining term 
we insert a factor 1 = dt'Cab{t,t')/dt'CabitTt'). We can then do another partial integral, and we could, in principle, 
repeat this process. All terms dX t — Tp share a common factor, namely the plateau value of Cabit,t'), 



C 



plateau 
ab 



(t) EE Cabit,t-Tp) 



(SMa{t)SAhit-Tp)) 



slow 



that also enters the relation 



1 

N 



{5Ma{t)5Mb{t)) 



fast 



N 



Cab{t,t)-Cf,^''''''{t) 



As a result, we derive from (2.34) our Ansatz (2.31) with a factor 

1 1 , dtfeit) Cabit,t') 



Te{t) fe{t) f^{t) dt'Cab{t,t') 



t'=t-Ti3 



+ 



This may be inverted, to yield 



Te(0 =Te(0+Te(0 



d\nCabit,t') 



dt' 



+ 



(2.37) 



(2.38) 



(2.39) 



(2.40) 



It is clear that the effective temperatures and Tg are not identical. However, in the models to be analyzed later 
on, we shall find that the difference is small. 

Notice that the ratio dt'C{t,t')/G{t,t') ~ Te{t') is allowed to depend on time t' . The situation with constant 
{— T/x, with X the break point of the Parisi function) is well known from mean field spin glasses, but we shall not find 
such a constant Tg in the models to be studied. In |15[| it is reviewed that in mean field spin glasses the fluctuation- 
dissipation parameter X{t,t')=TG{t,t')/dt'C{t,t') simplifies to X{t,t') = X{C(t,t'),t')= X{0,t') const. As our 
Te{t') will depend logarithmically on time, the i'-dependence of our X(0, t') = T/Te{t') cannot be neglected. We can 
therefore conclude that such a time-independence is an artifact of the mean field approximation. This supports our 
earlier conclusion that only at exponential time scales ~ exp(A^) the dynamics of the mean field spin glass is related 
to that of realistic systems [Q. In the numerical evaluation of the "fluctuation-dissipation ratio" T/T^ one should 
therefore keep in mind the realistic possibility of a slow time dependence of Tg. 



7. Time-scale arguments 

Consider a simple system that has only one type of processes {a processes) , which falls out of equilibrium at some 
low T. When it ages a time t at T = it will have achieved a state with effective temperature T^, that can be 
estimated by equating time with the equilibrium timescale. Let us define Te by 

t = Teq{T,) (2.41) 

We shall check in the models to be studied below that, to leading order in Int, it holds that = Te. (The first 
non-leading order turns out to be non-universal, since it already depends on numerical prefactors of T^q.) This equality 
also is found in cooling trajectories, when the system is well inside the glassy regime. It says that the system basically 
has forgotten its history, and ages on its own, without caring about the actual temperature. We feel that this is 
caused by the fact that dynamics in each new time-decade is basically independent of previous decade. 

In less trivial systems, for instance those having a Vogel-Tammann-Fulcher law, the timescale may have parameters 
that depend on the actual temperature, implying r^q — t{T,Tc,). We then expect that, to leading order, Tg follows 
by equating this expression with time t. 
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In many systems one finds a t' /t scaling in the aging regime of two-time quantities. There is a handwaving argument 
to explain that: 



(2.42) 

showing indeed the famihar t/t' scaling. In the models to be studied below we shall find logarithmic scaling corrections. 
They become strong at low T, and change the \fv~ft decay at T > to a t' jt decay at T = 0. So this argument might 
apply only to a subset of systems that fall within the scope of our approach. 



8. Results for simple systems that become glassy near T = 

In the remainder of this paper we shall consider two simple systems having only one type of processes [a processes), 
which fall out of equilibrium at some low T. Then the effective temperature Teify is expected to show up in the 
following deviations from the equilibrium situation: 

• matching the internal energy: U(t,H) — Ueq{Te{t), H) 

• matching the magnetization M{t,H) = Meq{Te{t), H) 

• from the configurational entropy via dU = T^dl — MdH 

• matching time with the equilibrium timescale: t ^ Teq{Te(t)) 

• via the fluctuation formula x""''* = Pe{t){5AP{t)) 

• from the fluctuation-dissipation relation: ^'^g*;* — Te{t')G{t,t') 

Even though these relations are not all independent, it is pretty clear that the whole glassy dynamics is strongly 
governed by one parameter: the effective temperature. 



3.. MONTE CARLO DYNAMICS OF UNCOUPLED HARMONIC OSCILLATORS 



Bonilla, Padilla and Ritort have recently considered an exactly solvable model with slow dynamics |2^ . It showed 
interesting, glassy behavior at low temperatures. In this section we present many details and further results for this 
model. This will also be a pcdogical step for the analysis of the spin model of section 5. 

After including an external field, the Hamiltonian reads 



^L=]^KY,xi-HY,^^ = \KM2-HM^ (3.1) 

i i 

where 



Mu = Nmu = V xt [k - 1, 2) (3.2) 



tt2 

E ^ Ne = n + N— (3.3) 
2K 



We introduce the shifted energy 



In a Monte Carlo step with parallel updates one replaces Xi ^ x'^ = Xi + ri/vN, defining = J^i^i + ri/y'N)'^. 
The thermal noise variables are independent Gaussian random variables with average zero and variance . For a 
parallel update of all Xi, this leads to the noise averaged transition probability from a state with (Afi,M2) to states 
with (M{, M^) = (Ml +yi,M2+ 2/2) 
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C nOO J ,2 \ 

n j ;7=pe" ^ 1 '^(^'^i - - 2^1 )^(A^2 - - y2) 

47rA V"^2 - TO? ''''^ 2A2 8A2(m2-TO2) ^ 
= P{yi,y2\mi,m2) (3.4) 

where we took the convention that probabihties involving extensive parameters are written with capitals, while those 
that involve intensive parameters are written in lower case. To derive this result, the delta-functions have been written 
in a plane-wave representation, and the limit N oo has been taken. Setting yi = y we introduce the variable x by 

K 2 
x = —y2-Hyi\ y2 = —{x + Hy) (3.5) 

A Monte Carlo move implies a change E' = E + x. The transition probability may be decomposed as 

P{y,y2\Mi,M2)dydy2 P{x,y\E, Mi)dxdy = p{x\e)p{y\x , e , mi)dxdy (3.6) 
with conditional probabilities 

p{x\e) = -=i= exp(-i^-^) (3.7) 

2 



p(y|a;, e, mi) = exp(--^^^^!^) (3.8) 



y 



having parameters 



xo = A!i^; = A^{K^m2 - 2HKmi + H^) = 2KA^e (3.9) 



-{KA^ - 2x){Km, - H) _ KA^ ~ 2x A^K\m2-ml) _ , K^ii 

~ 2{K^m2 - 2HKmi + H^) 4e " ~ K^m2 - 2HKmi + ~ ^ 2e ' ^ ' 

where we defined the deviations from equilibrium 

e ^^K{m2~ml + fil), = ^ - mi (3-11) 

z A 



We shall frequently encounter the energy scale 



(3.12) 



In one Monte Carlo step the probability of E evolves as [gOj 

P{E',t+^) = J dEP{E,t) J dxp{x\^)[W{f3x)S{E' ~ E~x) + {l~W{(]x))S{E' - E)] 

= P{E',t)+ J dEP{E,t) J dxp{x\^)WiPx)[6{E' - E - x) - 5{E' - E)] (3.13) 

here W{Px) = 1 for x < and W{Px) — e"*^^ for a; > is the Metropolis acceptance rate. The second term describes 
the rejected moves. Notice that the energy aspects arc independent of the field H, for the physical reason that H 
merely causes a shift of the equilibrium position, but not in E — (K/2) ~ H/K)"^ . This is due to the simplicity 

of the model. 

If one also keeps track of the magnetization, one has 

P{E', M\ t+j^) = P{E', M', t)+ J dEdMP{E, M, t) J dxdyP{x, y\E, M)W{f3x) 

X [S{E' - E- x)S{M' ~ M ~y)~ 5{E' - E)5{M' ~ M)\ (3.14) 
which, of course, is consistent with P{E, t) = J dMP{E, M, t). 
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1. Evolution of average observables 

We can now calculate evolution of physical observables. One derives from ( 3.13| ) 



{E{t+^)) - J AE'dM'E'P{E\M\t+^) = {E{t)) + J dEdxW{Px)PiE,t)xp{x\^) (3.15) 

where {E{t)) arises from the term without W{(3x). In the thermodynamic limit (i.e. for large N) P(E,t) will be 
sharply peaked around (E{t)), so one obtains a closed equation for the scaled average e{t) — (E{t))/N |po| 



dejt) 
dt 



dxW{Px)xp{x\e(t)) (3.16) 



This simplifying property is due to the lack of interaction between the oscillators. 
In the same way we proceed for the evolution of the magnetization 

1 r E E M 

{M{t + -)) = (Af(i)) + j dEdMdxdyW{(3x)P{E,M,t)yp{x\-)p{y\x, -, — ) (3.17) 

Here, and everywhere in the sequel, the y-integrals are Gaussian, and can be carried out analytically. This makes the 
problem with a field hardly more complicated than without. We obtain 

= J dxW{(3x)yop{x\e{t)) - ~{m{t) - |)/(t) (3.18) 

where 



fit) = - I dxW{Px)^p{x\e{t)) = I" dxW{Px)^-^p{x\e{t)) (3.19) 



2. Fluctuations 

The evolution of bilinear forms is a bit more involved. Let us consider the energy fluctuations. One has 

{E\t+^)) = (E-'{t)) + J dEdxWif3x)PiE,t)i2xE + x^)pix\^) (3.20) 



Using (3.15) this may be written as 



{EHt + ^)) {E\t))-2{E{t))mt + ^)) - {Em = J dEdxW{Px)PiE,t)pix\^){2xSE + x') (3.21) 
where SE = E — {E{t)). Expanding p{x\E/N) around E = {E{t)) one obtains for large N 

d{6E^)_ /'^^^(^^)((^2x|-+x2)p(^|£) (3.22) 



dt N J N de 

In the same way one derives for the evolution equation for fluctuations in M 

d (SAP) f , ,A5NP) d (SESM) So,,,,,, , , , 

^ - J dxdyW{px)i^j^2y— + ^ ^ ' 2y— + y')p{x\e)p{y\x, e, m) (3.23) 

The 2/-integral is Gaussian, and can be carried out. The result reads 

dt N=J d^^(/5-)(V^a;^2y„ + ^^-2yo + y^ + A>(x|e) (3.24) 

while for the cross-correlations 
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d {5E5M) r ^ ^^^^^ d (SESM)^ d . (SE^) d win., 

= / dxdyW{Px){ — T^x— + [x— + y—) + ——y— + xy)p(x\e)p{y\x, e,m) 



dt N J " ^' N dm N ' de "dm' N " de 

dxWmi^^ix^ + + ^-^'+xy,)pix\e) (3.25) 

A* OS dm N ae 



Recal ling that M — Mi and the definition (3.2) of M2, and adopting the definition of correlators to be given in eq. 
( 3.32 ), we may also cast these results in the form 

^C,t,{t,t) = J dxW{(3x){y^y, + Ay(^-^^ + ^ _^ (y^C,b(i, t) + y,C,a(t, i)) }p(a;|mi, ma) (3.26) 



where a,b — 1,2, Hi — H, H2 — —K/2, and, in the present model, p(a;|mi,m2) — p{x\£), with e defined in (3.11). 
Furthermore, 

_ x + Hiyo 

Vi =2/0 2/2 = (3-27) 

Since 5E ~ {K/2)5M2 — H6M1—— Hc^Mc, previous results are recovered from these relations. 



3. Correlation and response functions 



One can also consider the evolution of two-time quantities. The correlation and response functions for magnetization 
and energy are defined as 



Gmm (^1 ^ ) 


= ^ {5M{t)5M{t')) 


Gm7n ^ ) " 


Csm(t, t') 


= j^{5E{t)5M{t')) 


G^fji , ^ ) — 




= j^{5M{t)5E{t')) 




Cee{t,t') 


^ ±{SE{t)SE{t')) 





_ 1 S{M(t)) 
N 5H{t') 
. 1 5{E{t)) 
N 5H{t') 

. nt') 5{M{t)) 

N ST{t') 
Tit') 5{E{t)) 



(3.28) 
(3.29) 
(3.30) 
(3.31) 



N 5T{t') 

In eq. (|3.2[) we have introduced the macroscopic observables Mi = M and M2- The related correlators and responses 



are 



Cab{t,t')^j^{5Ma{t)5Mb{t')) 



_ 1 6{Ma{t)) 



(a, 6= 1,2) 



(3.32) 



They code the same information, but will be more useful at some stages. Cmm is just another notation for Cn, while 
Ce™ = (is:/2)C2i - HCii and C^e = {K^/4=)C22 - {HK/2){Ci2 + C21) + H^Cn. Similar relations hold for the G's: 



{K/2)G2i-HGii, G„ 



'iK/2)Gi2-HGii, and G,, = iA/K^)G22 + {HK/2)Gi2 - {HK/2)Gi2 + H^Gn, 



where we used that Gi2{t,t') = —G2i{t,t'). 

To derive the evolution of the correlations one considers 



(M(t + ^)M{t')) = J dE'dM'dEidMi P{E', M' , t + Ei, Mi,t')M' Mi (3.33) 
= {M{t)M{t'))+ J dEdMdEidMidxdyW{Px)P{E,M,t;Ei,Mi,t')P{x,y\E,M)yMi 

Subtracting {M{t + jj) — M{t)) {M{t')) and expanding P{x, y\E, AI), this yields the evolution equation for Cmm- The 
y-integral can again be performed. In similar ways one proceeds for the other correlation functions. One finally has: 



d f d d 

-g-^Cm^{t,t') ^ / dxW{Px)[C^mit,t')— + C,„,{t,t')—]yop{x\e) 



(3.34) 
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^Ce™(t,i')- f dxWif3x)xC,^{t,t')-^pix\e) (3.35) 



-C„,,(i,t')- / dxWipx)[C„,,it,t')— + Ceeit,t')—]yopix\e) 
ot J om oe 



(3.36) 



^Cee(i,t') = J dxW{Px)x-^^p{x\e)C,,{t,t') (3.37) 



Their equal-time values follow from the above fluctuation formulae. 
The equivalent formulation is 



d_ 

dt 



Cabit, t') = J2 Ccb{t, t')— / <:\xW{px)y^p{x\mi,m2) (3.38) 



4. Response functions 



The energy-energy response function Gee{t,t'), defined in (3.31), takes the form 

G,,{t,t')^^ J dEdE,dx^E[P{E,t\E, + xut' + ^) ^ P{E,t\E,,t' + j^)]^^^^p{x^\^)P{E,,t') (3.39) 

For our parallel Monte Carlo updates it holds that dt = 1/N. Both terms satisfy the same evolution equation, 
implying 

^G,,{t,t') = -g{t)G,,{t,t') (3.40) 
ot 

with 

git) = - J dx W{Px)x^p{x\e) (3.41) 

Since in the oscillator model the energy evolves independently of H, it is obvious that Gem = at all times. The rela- 
tion Gem — -~J2c ^cGci then implies G2i{t,t') = —{Hi/ H2)Gii{t,t'). In the spherical spin model, to be introduced 
later, this argument does not hold. 

From the evolution (3.18) of r7i(t) one gets immediately that Gmm and Gme satisfy 

^^Gmmit,t')^-f{t)Gmmit,t'); Gmeit, t' ) ^ - f {t)Gme{t, t') (3.42) 

The equivalent formulation is 

^Ga6(t, t') = J2 ^-bit, t')-£- j dxM^(/3a;)y„p(a;|mi, 7712) (3.43) 

The derivation of equal-time responses is a bit tedious. Let us take a = b = 1 and change the field from H to 
H + AH{k) at the time steps t + k/N, (fc = 1, • • • , tj). It holds that 

< Mi{t+'^) >-< Ah{t+-^) >= J dyidy2dMidM2WiP6Ek)yipiyi,y2\Mi,M2)piMi,M2,t+^) (3.44) 

where 

5Ek^^y2-{H + AHik))iM,+y,) + {H + AHik-l))Ah {k^l,---n) (3.45) 
Generalizing to all four case we have a response 
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, . 1 d<m,it + At)> 



Gab{t+,t)= lim — 



Tv^oo At dAHb 

,. < ma{t + ^) > -idem{AH{k) = 0) 
lim — 

N^oo dt J2k Ai/(fc) 

G':r{t+,t) + GT^\t+,t) (3.46) 



where we used that dt — 1/N and At = n/N . The main term arises from the j/^ terms in (3.45), 
G"r"(t+,t) = -/3 / dyidy2W'{px)VaVbP{vi,y2\m^,m2) 



= -pJdxW'{0x){y,y,+ (^-^^ Ay)p{x\m,,m2) (a, 6 =1,2) (3.47) 

where W'{Px) = — exp(— /3a;) for a; > and zero for x < 0. The contribution for switching on and off comes from the 
MaiAHaik) - AHa{k - 1)) terms, 

GZ''^'\t+,t) = j dy^dy2W'{Px)yaP{yi,y2\m,{t),m2{t))^ (3.48) 

does not depend on the precise switching procedure, but neither does it not vanish for adiabatic procedures. In our 
models this term will have contributions proportional to mi_2 ^ 1/t, so in the large-i domain it is much smaller than 
the terms of our interest. We shall neglected it from now on. We should point, however, out that the responses Gme 
and Gee do not involve such switching terms. This is related to the nature of the Monte-Carlo dynamics. 
For the responses with respect to an instantanous temperature pulse one has 



Gmeit+, t) = -(^J da; W'il3x)xyopix\e) (3.49) 
Gee{t+,t) =-[3 1 dxW'{Px)x^p{x\e) (3.50) 



The relation with Gab is exactly as for the C"s, see below ( |3.32[ ), as these two quantities donot suffer from switching 
effects. 

4.. GLASSY DYNAMICS OF THE OSCILLATOR MODEL 

We are considering a system with one mode. In view of the equilibrium relation Seq = T/2 we may introduce the 
effective temperature by 

Ut) EE 2e{t) (4.1) 

The dynamics of the oscillator model simplifies in the region T A — KA'^ /8. Technically this occurs since for 
£ ^ A OT Te A we can approximate in the expression 

the last, Gaussian factor by 1 — x^ /WATe, leaving only exponential integrals. We shall investigate the dynamics in 
this region, and look for relations satisfied by observables. 

1. Equilibrium regime 

In equilibrium one has e ~ T/2, rn = H/K. Then there holds detailed balance 

W{Px)p{x\e) = W{-(3x)p{-x\e) (4.3) 



17 



assuring that e = in eq. (3.16). For eq. (3.22) this implies 



in accordance with the equihbrium relation AU/dT = l3^(SH'^). The relation {K2§ amounts to 



(4.4) 



(4.5) 



This is also expected, since only the diagonal i — j terms in (SM^) = ^^ASxiSxj) contribute at equilibrium, showing 



that K {SM"^) /2 will indeed reduce to NT/2, the equilibrium value of E. Finally, eq. (3.25) tells that the cross- 
correlation {5E5M) vanishes. 



The evolution equation (3.16) for the energy can be expressed in the notation of ref. 



Ae(t) 
dt 



-(2£(i)-|)/(t) + 2A erfc(a(t)) 



where A = KA'^/8, a{t) = y/A/2e{t), f{t) was defined in eq. ( p9| ), and 



erfc(Q!) 



S^^-2^ + 4^) 



(4.6) 

(4.7) 
(4.8) 



We can look at relaxation close to equilibrium, where a^^q = ^/ (3A. We set e — T/2 + fe. Eq ( ^.16| ) becomes to 
linear order 



^ = -5e \8(iA{2j3A + l)erfc(v/^) - 16/3AW ^e"^^ ) = 



dt 



Se 



(4.9) 



When T ^ A the equilibrium timescale becomes, due to the expansion (4.7) 



7r(3A 



64 



(4.10) 



The latt er als o follows from ( 3.16 ) by performing the a;-integral, after neglecting the exp(— a;^/2A2;) factor of p{x\£) 
From (3.18) it is clear that the magnetization relaxes to its equilibrium value rrieq = H/ K at timescale 



j{m) 
'eq 



2T 



/(oo) A 



(4.11) 



The important fact for us is that both time scales have an Arrhenius behavior ^ exp(j4/T) at low T. This implies 
that the oscillators, subject to parallel Monte Carlo dynamics, can easily fall out of equilibrium at low enough T, and 
thus exhibit interesting glassy behavior. 



2. Cooling procedures and the glassy transition 



Eq. ( |3.16D simplifies in the regime <ti A. Indeed, as ~ ^ A^; = SAT^, we can neglect the Gaussian factor 
exp(-x^/2Aj,) of p{x\e) in (O). We thus obtain 



2T? 



:(1 



rp2 



V^fATe' (2re-r)2 



(4.12) 



The condition Tg > T/2 is typically satisfied, since T^ > T in cooling and aging experiments, and Te ^ T in heating 
the glass. 



Using (4.10) we can write this for Tf, close to T in the universal form 
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T-T, 



' eq 



(4.13) 



We can now introduce the inverse function T^g^{t); in our case ( 4.10 ) it reads to leading order T^^{t) — A/h\t. Let 



us then consider a non-linear cooling process of the form llC 



■ eq 



Tit) = (1 - R)Tg + RT-\t) « (1 - R)T, 



(4.14) 



It involves two parameters: the glassy transition temperature Tg and the dimensionless cooling speed R. A non-linear 
cooling experiment of this form could be performed in any system with a quickly diverging equilibrium timescale. 



We first show that a glassy transition occurs around time scale tg — exp/S^A, where one has T{t) 
timescale during which the system basically remains at temperature T is 



'^cooli 



m 
\m\ 



Its ratio to the equilibrium timescale is 



'^cooling im) ft 



' eq 



im) 



t 



R-l 



T„. The 



(4.15) 



(4.16) 



We can discriminate three cases: a) When R > 1, then for t <^tg there is equilibrium at the instantaneous temperature 
T{t). For t> tg the instantanous equilibration time r^q is larger the cooling timescale Tcooiing, and the system becomes 
glassy, b) For < i? < 1 this process describes cooling in a glassy state so slowly, that equilibrium is reached around 
time tg. c) Finally, for i? < it describes heating in the glassy state, and equilibrium is reached around time tg. 



Eq. ( [4.14D imphes that 



R 



T',g{Tg-{Tg-T)/R) 



This allows us to combine (4.13) and (4.17) into the time-independent form 



dT 



T-T, Tlg{Tg-{Tg~T)/R) 



H 



R 



' eq 



(Te) 



(4.17) 



(4.18) 



For T ~Tg ^ Tg /A one has the equilibrium value Tg — T, with expo nentia lly small corrections. Well below the glass 
transition T-Tg < -T^/A one has = Tg-{Tg~T)/ R. Due to eq. ( [l.lO|) this imphes that TeajT^) exp /3eA - so 
Te K, A/ Int. As we shall see in next section, this is the same behavior as occurs for aging at T = 0. We may conclude 
that the system then basically has forgotten its cooling history, and just ages as at any low enough temperature. 
Similar behavior was found by Godreche and Luck in the backgammon model . 
We may go to dimensionless variables by putting 



T^Tg + -j-x; 



T, = T,. 



A ' 



(4.19) 



and obtain 



dy ^ y-x ^y_^/p. 
dx R 



(4.20) 



This equation is probably universal. Indeed, it is a small excercise to check that the very same equation follows from 
(4.18) when r^g has a Vogel-Tammann-Fulcher-type law T^g ~ exp(A^(T — Tq) '), and a glass transition occurs in a 
narrow range around some Tg with Tg — Tg <C A. 



Equation (4.20) can be solved analytically for large negative and large positive x. Let us introduce 



(4.21) 



For large negative w one sets 
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y = x-\iiz{vy, v = e'" =e-^^-^^''/^ (4.22) 

to obtain 

{R - l)v^z'{v) + hiz{v) + Rvz{v) = (4.23) 

By series expansion one finds 

z^l~Rv + —(5R- 2W (29R^ - 27R + 6)v^ + — (llSli?^ - 18127?^ + gooi? - 144)z;^ 

2 3 24 

R 

- -(1529R^ - 3345i?2 + 2690i?2 - MOR + 120)w^ + • • • (4.24) 
5 

This implies for tlie specific heat factor an exponential approach to equilibrium 
dy lnz(z;) 



dT 



H Ax Rvz{v) 

= 1 + {R-1) (-v + 2{2R-l)v^ - h5R-2){3R-2)v^ + ^{iR - 3){29R^ - 27R + 6)v'^ 



(4.25) 



When i? > 1 or i? < this applies for large positive x. If < i? < 1 it applies for large negative x. 
For large positive w one sets 

'^- = -7W^ (4.26) 

and 

y^-—l— + lns-\nu{s) (4.27) 
(R — l)s 

yielding 

u{s) = 1 + sins - slnu(s) - (i? - l)su(s) + (i? - l)s'^u'{s) (4.28) 
By iteration one finds an expansion in powers of s and A = In s 

u{s) = 1 + {A - R+l)s + {2R-2- A)s^ + {-8R + y + 3A + ^A^ - 2Ai?^ + ^R^)s^ 

+ i-^R" + ^fi' + 45i? + 18Ai? + 2A2i? - 5Ai?2 - ^A^ - Ja^ - 14A - ^)s^ (4.29) 
2 3 3 2 6 

This implies for the specific heat factor an algebraic approach to the value 1/R, with logarithms in the subleading 

terms 

dTe _l + sA - slnM(s) 
'dT Ru{s) 

= 4 + ^^--^ [s + {R- A- 2) s^ + {R^~2RA~ 7 R + A^ + 7 + 5 A) 
R R 

+ (2i?3 - 35i?^ - 6i?^A + 46i?A + 92i? + GA^R ~ 52A - ITA^ - 61 - 2A''')y] + • • • (4.30) 

When i? > 1 or i? < this applies for large negative x. If < i? < 1 it applies for large positive x. It is trivial to see 
that both (4.25) and (4.30) go the correct value dTe/dT = 1 in the equilibrium limit R ^ 1. In that limit the system 
will remain in equilibrium, because the cooling procedure is very slow. 

In figures 4.1-4.4 we present the universal lineshapes for the specific heat factor dT^/dT ~ 2c for several values of 
R. They exhibit the features known from experiments. 
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A(T-Tg) / 

FIG. 4.1. Specific heat factor dT^/dT as function of reduced temperature in non-linear cooling experiments with different 
speed parameter R. The asymptotic values are 1 to the right and 1/R to the left. 
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FIG. 4.2. Specific heat factor dT^j&T as function of reduced temperature in non-linear heating experiments with different 
speed parameter R 
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-10 -5 5 10 

A(T-Tg)/Tg2 

FIG. 4.3. Specific lieat factor dT^/dT as function of reduced temperature in a non-linear cooling ex perim ent with R = 2 and 
in a non-linear heating experiment with R — —2. Dashed lines are asymptotes from eqs. ( 1.25 ) and ( 4.3C ). 
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T/ A = 0.045, R = -3 
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FIG. 4.4. Specific heat factor dT^/dT as function of reduced temperature in a non-linear cooling experiment with R = ?>, 
and for a non-linear heating experiment with R = —2>. The "universal" curve is from previous theory. In the "real model" the 
original differentio-integral equations have been used, and R has changed from ii = 3 to J? = —3 after reaching T = 0.045A. 



The analysis of this section thus shows that cooling in systems with an Arrhenius law leads to glassy behavior quite 
similar to that expected for realistic glasses. 



3. Aging in the glassy regime 



Suppose we quench the system at time t — Q from an equilibrium state at temperature Tquench ^ A to a lower 
temperature T. As opposed to previous section, we now assume that after the quench the system is far from equilibrium 



(viz. Tquench - T > T^^gnc/iM)- Then eq. ( [4.13D may be written as 



where we use the short hand 



(l-r2)V^ 



T 



-At 



2T,-T 



(4.31) 



(4.32) 



22 



Its integral is 



where to follows by inserting at t = the value Tg = Tquench- This result may be written as 

/3e A - i In /3e A = In ^ ; tq xt/^ (4.34) 
2 To 2(1 - r^j 

For large t one has to leading order — A/\nt, while the initial condition gives a small correction of order to/t ^ 
exp(— /3eA); it may thus be neglected for t > lOip- This says the initial condition is washed out, and is the basis for 
our interpretation that each decade is practically independent of the previous one. Likewise, the effect of a finite T 
is very small, and to leading order one could set T = 0. This says that in the glassy regime the energy essentially 
evolves as if the system had been quenched to T = 0. Only near the return to equilibrium the T-dependent factor in 
34h brings a vanishing argument in the logarithm, from which non-trivial behavior results, as discussed in section 



To leading order one may invert eq. ( 4.34 ), to obtain 

To 2 To 

In practice this need not be a good approximation since 1/lnt is usually not very small. For our purposes (leading 
order expansion in powers of Tg) this is equivalent to A/ {\nt + i Inlnf), and actually even to A/ Int. It is a simple 
excercise to check that one has 

r^nTe)^il-J^^-Jy)^-t (4.36) 

proving our general assertion that in the aging regime the effective temperature also follows by equating the equilibrium 
time scale to t. 

In the aging regime eq. (3.19) becomes 



/ = 2^e-P'^{l + r)(l - ^) = 7^ = (4.37) 



Te.r'^ 1 _ 2Te - T 

~2A' ^ {l-r)t ^ 2{Te - T)t 

Let us define its integral as 



Using dln/i/de = //e we obtain 



ft 

h{t) = const X exp / dt'f{t') (4.38) 



h(t)-(i m.jM)!^!£!:^ .439) 

Ht) - (1 - — ) _ (4.39) 



For T > this behaves as ^/i, with logarithmic corrections. At T = the resulting asymptotic scaling h{t) ^ t(\nt)^ 
differs from the result h ^ t\/lnt reported by |20[ . However, their figure 2 already show a deviation between the 
data and their asymptotic formula, that becomes increasingly worse in the asymptotic limit. We were informed by E. 
Hennes that the present expressions ( 4.39D , (4.34) give for t' > 100 an almost perfect agreement with the numerical 
solution of the integro-differential equations 1 39| . 
For later use, we mention the results 



g = 2e-^^^J^{l -r' + ^r\l + r)(l + 3r)) (4.40) 



and 



h{t) = const X exp /* At'g[t') = [P^Af'^ ^ ^^^'^^^ e^'^ (4.41) 

Jo t — Pe 
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Due to eq. (3.18), the magnetization relaxes as 



In the regime of large times and small ~ l/\nt, the deviation ■m{t) from H/K is exponentially small in Tg. As 
compared to the powerlaw that occurs in the energy (recall that e = Te/2), this can be neglected. This says that the 
magnetization quickly goes to its quasi-stationary value. 

4. Correlations 

Bonilla et al. considered the on-site correlation function {xi(t)xi{t')) At non-zero field it would become 

{5xi{i)6xi{t')) ^ with 5xi — Xi— <Xi> (t). However, the correlation function related to thermodynamics is the global 



correlator Cmmit,t') = j{Sxi{t)Sxj{t')), defined in (B.28). Its equal time value is found from eq. (3.24). We shall 
now study its two-time structure. 
Let us introduce 

f dxW(f3x)xp(x\e) T^-T , , , , 



where the explicit result holds for and T much smaller than A. After dividing by ( ^.l(j ) we may write ( 3.24 ) as 



2<x> ^ ^\ 2<x> -KA'') + (4.44) 

where we used that < SESM > is exponentially small in T^. The dominant behavior follows by neglecting the < x >~ Te 
terms. However, they can be fully taken into account, as the solution to this equation reads 

C„m(i,i)-^ (4.45) 

Corrections are exponentially small in Te- This result even holds when T, occurring in < x >, depends on time. It 
also allows a simple check: The result K {SNP) /2k,E = NTe/2 is in accordance with the expectation of Bonilla et al 
that off-diagonal terms (SxiSxj) with i j are subleading. 



Solving eq. (3.34) for t' we can neglect the C^m term, since it is exponentially small in T^. This yields 



Te(t') h{t') 

Cmn.{t,t) = ^^ (4.46) 



where h(t) was defined in eq. 



a. Stretched exponential fitting and arguments against doing that 

In the study of glasses, where often at best two orders of magnitude of C can be determined, it is commonly assumed 
that there occurs a stretched exponential decay, 

Crnmit, t') = a{t') exp(~(i/r)^) (4.47) 

It is often stated that such stretched exponential decay is one of the basic properties of the glassy state. 
In our case we would need that 

/i(i) = exp(-d-|- (i/r)^) (4.48) 

for some set of parameters d, r, 7, or, equivalently, that ln(c? -f ln/i(i)) is linear in Int with slope 7. In view of the 
exact expression, this is certainly not an exact description. Let us, however, look at the plot for d = in figure 4.5. 
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FIG. 4.5. The decay of the correlation function can be described by a stretched exponential in a time-window where the 
data for Inln/i are linear in Int. The stretching exponent then equals the slope in this figure, and will depend on T and on 
the chosen time-window. From left to right: T/A — 0.1, 0.075, 0.05, 0.025, 0. The bending of in the lines on the left (having a 
relatively large T) indicates that equilibrium is approached at the considered timescale. 

It is seen that lnln/i(<) = 7lnt — 7lnT can be a reasonable approximation in a not-too-wide large-time window 
tmin < t < tmax- In agreement with usual findings, the effective exponent 7 will decrease with T, and be bounded by 
the finite T — Q value. Notice, however, that it will also depend on the time window where the fit is made. 

In the stretched exponential fitting procedure there is one more adjustable parameter, namely the overall prefactor 
exp(— d). In figure 4.6 we take T = 0.0025 A and give plots of ln(c? -)- In ft,) versus Int for various d. In intervals where 



this curve is flat, h is well described by a stretched exponential (4.48) 



T=0.0025A 

ln(d+ln(h)) 
5 d=100 




10 15 20 25 30 
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FIG. 4.6. Log-log plots of the function d + lnh{t) at T = 0.0025 A for various d. h is well described by a stretched exponential 
h{t) ~ exp(— ci+ {t/T)'') in an interval where one of the plotted lines is straight. Then the slope yields 7 and the offset — 7lnr. 

This information can be used to obtain C{t,t') ^ h{t')/h{t). To exagerate what happens, we take a very simple 
linear fit to the data of figure 4.6: we consider the interval 10 < Int < 20 and make, for a given value of d, a linear 
interpolation through the data points at Int = 10 and Int = 20. Taking t' = exp(lO) we plot in figure 4.7 t he fits to 



\og{hlt') / h{t)) for the cases d = \Q, d = 100, and compare with the exact result h{t')/h(t) from eq. ( 4.39 ). By our 
construction, the results agree at \nt — 10 and 20. It is seen that in all cases the fits are reasonable in regard of the 
scale presented in the figure, and that increasing d improves the overall fit. 




FIG. 4.7. Log-log plot of stretched exponential fits to the ratio h{t')/h{t) at fixed t' — exp(lO) and T = 0.0025 A. Upper 
line: by taking a linear interpolation between the points \nt — 10 and Int = 20 for the case d = 10 of previous figure. Middle 
line: the same for d — 100. Lower line: exact result. 

The free fitting parameter d is not present in reality. It occurs if one overlooks that h{t)/h{t') should be fitted as 
function of two parameters, namely both t and t' . Indeed, at one given t' one is free to choose d; however, emposing 
the asymmetry t ^ t' would bring d {t' /t)'' . This reduced freedom would decrease the overall accuracy of the fit. 
In practice this is typically not done, partly because of lack of data curves C{t,t') at enough t' values. Nevertheless, 
whenever fitting of stretched exponentials is attempted, we stress to make a two-parameter fit of C{t,t'). 

In our opinion the present fitting procedures merely say that stretched exponential fits can too often be made, 
without yielding too much insight. The standard statement that stretched exponential decay is one of the basic 
properties of the glassy state should, in our view, be taken with a big grain of salt. What really happens is a slow 
decay, of which too few orders of magnitude are known to draw firm conclusions on their analytical form. It seems 
needless to say that in experiments the tails of the correlations have large relative errors, which, in our opinion, make 
the problem quite insensitive to stretched exponential or many other fitting procedures. 

We feel that the situation even becomes worse if such fitting is applied for showing the presence of critical behavior 
near a supposed critical temperature in the glass, as is sometimes done. 



5. Fluctuation-dissipation relations 



We now consider aging dynamics at fixed T <ti A. Neglectini 
time value Gmm{t,t) — f{t)/K. This result is exact for i7 = [| 



^all j/o contributions, we find from ( 3.46 ) the equal 
!0[ . Its two-time form follows from ( 3.42 ) 
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Gmm (^1 ^ ) 



fit') hit') 



K hit) 

We can now consider the fluctuation-dissipation relation. We define the effective temperature Tg by 

dt'Cmmjt, t') _ - , 

By direct evaluation of the left hand side, we find from ( 4.46| ), using (4.12) and (4.37), 

Teit) 



TS) = T,it) 



fit) 



^ A 2Te-T 



(4.49) 



(4.50) 



(4.51) 



Due to eq. ( |3.34 ) this agrees with the general result ( ^.40| ), which assures consistency with single-time expressions. 
As the second term is smaller by a factor T^/A with respect to Te, we thus see that to leading order the same effective 
temperature occurs as in the energy and the time-scale. (The same happens when the local correlator is considered 
pO|). Notice, however, that the leading correction is non-universal, as it depends on the model parameter A. When 
equilibrium is approached, — > T, and the correct limit T,, — T is reproduced. 

Let us now look at the energy fluctuations. To leading order we may neglect d{SE'^)/dt in eq. ( p2| ). This brings 



N 



<SE^>^C\,it,t) 



-T^<x^> Tj^it) I 



4:A<X> 



A 1 -r2 



(4.52) 



This result shows various points. First, even at r = T = it is one order of magnitude smaller than what one would 
anticipate from the equilibrium expression C^e = Only in the glass transition region T^ — T^ (1 — r)T ~ T'^ /A 

the equilibrium scaling C^e ^ is recovered. We conclude that there is a complicated, non-universal T-dependence 
in the whole aging regime T < Tg. The possibility of a model-independent generalization outside equilibrium of the 
relation dJ7/dr = 0^{5E'^) will be discussed in the Discussion. 
For different times we find from (3.37), (4.41) 



Ceeit.t')^ 



T^it') l + r^it') hit') 



A l-r2(t') hit) 
Likewise, the Greens function follows from ( 3.50| ) and ( 3.4C| ) as 

G,,(t,0 = i^^'^ ^2r,V(i + ^)e-/5.A 



hit) 



if) 



To leading order this yields the fluctuation-dissipation relation 

dt'Ceeit,t') _ rpiee)/,/-. 
Geeit,t') 

where the t-dependence again has dropped out. The quantity 

re(4T2-6Ter + 3r2) 



(4.53) 



(4.54) 



(4.55) 



r(ee)(,. Ut)il + rHt)) 

' ^' r^t)il+rit)) 



2^2 



(4.56) 

T. Unlike Tp itself, it has no obvious 



is also an effective temperature that also has the correct limit when Tg 
model-independent interpretation. We feel that this is due the fact that it relates to subleading quantities. 

For the fluctuations ( 4.45 ) in M the non-universal terms are exponentially small in Te- This shows that for the 
fluctuations in M = Mi a quasi- universal behavior takes place. As Cem is negligible, the M1M2 cross-fluctuations are 
also simple. 



Cl2it,t') 



2HTeit')hit') 
KV^) 



(4.57) 



plus exponentially small corrections in Tg, or powerlaw in 1/t. The M2 correlations and responses have the form 

4/^2 4 att2 a 

C22{t,t') = — Cn(i,i') + —C,,it,t'); G22it,t') - — Gn(t,t') + —G,,it,t') (4.58) 
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In both expressions the first term is two orders of magnitude in Te larger. On top of that, the second term decays 
faster (~ t' /t versus \fV~ft) whenever T is non-zero. It thus holds that in all four cases 

K&=1,2) (4.59) 

This simple result suggests that in general the fields could also stand for a chemical potential, a pressure, or a quenched 
randomly directed forcing strength. 



6. Non-equilibrium thermodynamics 

We now wish to view previous results in the thermodynamic framework of Section II. 

Since there is only one type of processes, that are by definition the slow ones, the entropy of equilibrium processes 
Sep vanishes. For such cases the configurational entropy can be derived simply. It is defined by the degeneracy of 
states with energy [/, and given by the micro-canonical partition sum 



= J Dx 5{n{x) -U)= j ^ ^ y e-''«^"^ (4.60) 

where Dx — YliAxi is the integration measure. By the saddle point method one obtains 

I = max /?([/ - F,q{f)) = 5eg(Te) (4.61) 
f 

where we used that for U = Ueq{Te) the minimum is assumed at T = T^. This result holds generally in simple systems 
with only one timescale that diverges near T = 0. Here we have 

T 

X=5e,(re) = y(l0g-^ + l) (4.62) 



Since dU — NdTe/2 — {NH/K)dH, it is now clear that the formulation (2^) of the first law is satisfied in the present 
non-equilibrium state. As S'op = 0, the free energy reads 

F = U -TJ (4.63) 



and it satisfies the relations (2.11) 



As mi = H/K is temperature- independent, the modified Maxwell relation redu ces to the standard one: in eq. 
( ^.15 ) the terms proportional to T vanish, and the ot her t erms follow already from ( |2.9| ) with S'op = 0. Neither is it 



interesting to investigate the first Ehrenfest relation (2.21): It holds trivially, as one has a = 0, x ~ ^/K, implying 
Aa = Ax = 0. Notice, however, the present results already require that second Ehrenfest relation ( ^.22| ) is modified 



outside equilibrium ||6|. Indeed, from section ^ we have AC = N{R - 1)/{2R) ^ 0, while Aa = 0. Equation ( |2.24| ) 
is nevertheless satisfied, since DT^/dT = and dJ/dT = dl/dT = N/{2Tg) 



The fluctuation formula ( ^.31 ) are also satisfied. To show this explicitly, let us take a — b = 1. Since there are no 



fast processes, the first term vanishes. The same holds for the third term, since m = H/K leads to dm/dTe = 0. Due 



to eq. ( 4.45 ) the second term equals {NTe/ K) / {NTe) = l/K, which is the desired result. We can also check it b y 
integrating the up the instantaneous field pulses, as was done more generally in the argument starting with eq. ( ^.35 ) 



The same conclusions hold for the other three cases. We have already mentioned that the fluctuation-dissipation 
relation ( ^.34 ) is satisfied with Tg given in eq. (4.51), and that the apparent specific heat Ch = N{dTe/dT)/2 has 



no simple connection with the energy fiuctuations ( 4.52| ) 



In all situations considered we have seen that Tg ^ A/ hit is to leading order in agreement with the timescale 
relation Teq{Te) ^ t. We have also seen that correlation function has the scaling h{t')/h{t), with, at finite T, 
h(t) ~ times a function of Int. In the T limit they become so strong that they replace the h{t) ^ \/i scaling 
by /i = t X function(lnt) 

In conclusion, the proposed picture applies to the harmonic oscillator model, be it that a few aspects are quite 
trivial. In next section we shall consider a model of spherical spins, which has a richer behavior when changing field, 
and in regard to the Ehrenfest relations. 



28 



5.. MONTE CARLO DYNAMICS OF FREE SPHERICAL SPINS IN A RANDOM FIELD 



Previous model had the drawback that the effect of a field was rather trivial. It was therefore of no great interest 
to check the first Ehrenfest relation: it is satisfied in a trivial way, having Aa = Ax = 0. 

We have therefore considered a closely related model, containing free spherical spins in the presence of a random 
external field, which does not share these drawbacks . Also in this very simple model Monte Carlo dynamics can 
be solved exactly and leads to glassy behavior. In fact, the dynamics just maps to leading order onto the one of the 
oscillator model of previous section. 

The Hamiltonian contains two parts, a "self- interaction" term involving fields F^, and a coupling to an external 
field H 

N N 

n=-J2'^A-Hj2s., (5.1) 

i=l 2=1 

The model is solvable for any set of quenched random fields Ti that have average zero and variance T^. To simplify 
the discussion, we make the additional (but technically unnneeded) assumption that Tt — ±T, implying that at each 
spin position there is a quenched random unit vector Ti/T, along which the spins wish to point for large pinning 
field r. This limitation allows the exact gauge transformation Si — > TiSi/T, which interchanges the role of H and F. 
Without the additional assumption this interchange would also exist; this is due to the spherical nature of the spins. 

In terms of the "staggered" magnetization Ms = (1/F) J^i ^iSi one simply has H = — FM^ — HM. When defining 
Hi ^ H, H2 = F, Ml = M, M2 = Ms, we may also write this as 

2 

n = -Y, HcM, (5.2) 

c=l 

The spins are spherical, which means that they can take all real values compatible with 

Y^S^^^N (5.3) 

i 

In equilibrium the system has a free energy 

— f = — log 3n - ^ (5.4) 

where 



K = v/F2 + m (5.5) 
and with chemical potential fixed by optimization, implying 



^, = ^jK^ + }-T^ + hr (5.6) 



This yields for the internal energy, the magnetization, and for the entropy 



1 H H HT HT^ , , 

— Me„ = — w H (5.8) 

iV ^ ^ K 2K^ 8K^ ^ ' 



1 I T 1 I T 1 

— 5e„ - -ln- + - w -In— + - (5.9) 



The approximations hold for low T. 



29 



1. Monte Carlo dynamics 



As for the oscillators, one makes parallel Monte Carlo moves Si S[ — Si + ri/\/N, with the independently 
drawn from a Gaussian with average zero and variance A^. Next one makes a global rescaling of the length of the 
spins, to reinforce the spherical constraint. This leads to the final update per time step 



Sj — Si 



ri_ 
N 



(5.10) 



This conserves the constraint (5.3). It implies for the change in the energy and in the total magnetization M = Si 



(r, + H)r 



N 



-n- 



NVN 



n 



27V2 



}; M' - M = ~ 



M 



(5.11) 



Introducing the new variables 



M 

TO=— ; fii = meq{e) - m 



H H 



(5.12) 



which are small near equilibrium, this leads to moves Ti' = Ti. + x, M' = M + y with Gaussian transition probabilities 
of the type (3.6), (3.7), having parameters 



(5.13) 



2/0 



H 



A^if 2 _ 2Kx + 2xs 
2e{2K-e) ' 



e{2K-e)' 



(5.14) 



where 



K = + m 



(5.15) 



In particular at small e the model of spherical spins in a random external field leads to a problem very similar to that 
of the uncoupled, identical oscillators in a steady field. The previous general formulae for variances, correlation- and 
response functions remain valid here. 



2. Glassy transition and the Ehrenfest relations 



The evolution for e again satisfies eq. ( 3.16 ), with the new expressions for xq and A^:. We can again introduce Te 
by equating U = —K + e — Ueq{Te), so that 



T.= (^^^.2e(l 
K-e ^ 



2K 



); 



e = K 



2 



8K 



We also define 



A = 



8 



A2 

B=—{K^e) 



(5.16) 



(5.17) 



p{x\e) takes the form (4.2) with A ^ B. For Tg <^ B one can again approximate it by an exponential. This yields 
the equilibrium timescale 



'eg ~ ' = ^— 1 ' ^ ■ ^ 
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(5.18) 



For a non- linear cooling process of the form (4.14) the results of section 4. 2 apply immediately. We consider cooling 
sequences with i? > 1, where the system goes from a paramagnet to a glassy state in a region around some Tg <^ B. 
Below the glassy transition one has an apparent specific heat 
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c = a 

This is of the general form (2.8), with background Ci = 
holds that 



dT 



H 



(5.19) 



0, since there are no fast processes in the present model. It 



C2 = 



N 



4^X2 + T2/4 



4K 



(5.20) 



The same universal lineshapes of the oscillator model thus occur here, with non-trivial prefactor C2 . Between the right 
and left sides of the glassy transition region there is a difference AC = A^C2(1 — 1/^)- 

Furtheron it will become clear that m{t) — meq{Te) remains zero upon cooling at fixed field H. Therefore one has 



m = meq(Te) 



This yields a magnetizability 



H 
K 

H 



H _ H 



if 2 



4 2 ' 



(1- 



dT, 



and, since = r2 + ^ a susceptibility x = Xii of the form (2.28|), with 



where w = yjK"^ + /4, and 



.fluct 



conf 



r2 + r2/4 + wTj2 



H 



A"3 



H 



X 



w{2w + Te) dH 2K 



(p2 _ ff2)j.^ 

2K'^ 



2*- 2K' dH 



(5.21) 



(5.22) 



(5.23) 



(5.24) 



Around the glassy transition there occur smeared discontinuities in the apparent specific heat, magnetizability and 
susceptibility 

dT 



Ac= C2(l 



.) 



]pC2(l- 

H dT, 



H 

dT, 
dT 



H 



These results and eq. ( 2.25| ) allow us to verify the fir st E hrenfest relation ( p. 21 ). 
From the identity T = S,q{T,) and the expression (5^) we obtain 



dl 
dT, 



T.H 



9i 

T,' 



dl 
'dH 



HC2 
if2 



(5.25) 
(5.26) 
(5.27) 

(5.28) 



We can now consider the modified second Ehrenfest relation ( 2.24 ). Due to eqs. ( |3.25 ) and ( ^.26 ) it takes the form 



C2_ 

T, 



(1 



dT, 
dT 



H 



He, 

K2 



dT, 
dT 



dH^ 
h' dT 



(1 



dT, 
dT 



.)( 



H 



C2 dT, 
Tg dT 



Hc2dH„ , C2 dT, 



H dT 



Tg dH 



dHg 

T dT 



(5.29) 



After dividing out the common factor 1 — dT,/dT and eliminating the remaining dT,/dT by use of the relation ( ^.25 ), 
it is seen that terms with and without dHg/dT cancel separately. This implies that the modified second Ehrenfest 
relation is satisfied for any value of that parameter, as was to be expected. 

Using ( ^.27 ) and the relation Tg « A/ Intg, the Prigoginc-Defay ratio can now be expressed as 



n 



AC dT, 



9 



TNAa dH 



dT„ 



HTg dH 



= 1 - 



if 2 dlnt„ 



ii In i„ dii 



(5.30) 



Contrary to what was long believed, the condition 11 < 1 is easily met. Indeed, in case that R is fixed, we may 
still choose the glassy transition line Tg{H), or, equivalently the glassy transition timescale tg{H). Values 11 < 1 
thus occur when dtg/dH > 0, so in half of the sets of smoothly related cooling sequences. This analysis confirms 
our general argument that the Prigogine-Defay ratio can take any value between zero and infinity, and perhaps even 
negative values. 
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3. Aging regime and its thermodynamics 



For temperatures in the aging regime we have, very analoguous to (4.31), 



dt [K 



(5.31) 



where we again use the short hand r = T/{2Te — T). To leading order this relation even reduces to eq. (4.31), the 



only change being exp(/3e^) exp(/3eS) ~ exp(/3eA + A/2K). Therefore to leading order again follows from eq. 



( [4.34D , with A + ATe/2K. 

Let us stress that we donot consider the regime A ^ 1, where a non-universal regime 1 <C Int <g; would occur. 
This is the subject of a recent work on a related model with fast and slow processes, in which a Kauzmann transition 
occurs . 

In next subsection it is made clear that jii = meq — m is exponentially small in Tg. Therefore the magnetization 
very closely follows its quasi-equilibrium value set by Tg. 
We can now check the thermodynamics. It holds that 

U = N{-K+^-^) (5.32) 
M = N{--^ + ^) (5.33) 

X^N(-\n^ + ---!^) (5.34) 
4 K 2 AK' ^ ' 

The first law ( |2.9| ) with Sep = reduces to dC/ = Tgdl — MdH. At constant H it holds because one has replaced 
T — !■ Te in energy and entropy. Then one can take the difference between aging experiments at two nearby iJ's. The 
relation remains valid because M is essentially equal to its quasi-equilibrium value M^qiT^,, H). 



In the modified Maxwell relation (2.15) the terms without explicit factor T cancel because of the first law ( |2.9| ) with 
Sep = 0. The terms with explicit factor T cancel separately, because of the quasi-equilibrium relation &I j dil\T = 



—dM I dTf\T,H that follows from (2.11). Such a pairwise cancellation would, of course, not occur in less simple models. 
Indeed, in a glassy model with directed polymers the modified Maxwell is also satisfied, but in a non-trivial 
manner |^. 

4. Changing the external field and the need for an effective field 

If also field H is changed in time, then we have to be more careful. The Monte Carlo dynamics now leads to 

ma= \ dxW{l3x)yj>{x\e) (5.35) 



where, again, m\ = m, m2 = vria and i, = x — Hm. We shall be interested in cases with logaritmically slow H, 
implying H l/t, where x ~ x. 

From the definitions ( 5.12| ) we derive for the evolution of e and /ii 



e = H{— -m)+ dx W{Px)xp{x\e) = H{^e + ^ii) + / da; W{f3x)xp{x\e) (5.36) 



Ai = -H{-^ + e— ^ + (- - m) — ) - /Ml -^H—~ M (5.37) 



So far we considered cooling at constant field. Then eq. ( 4.39| ) says that fJ-i{t) — fj,i{to)h{tQ)/h(t) decays as a 



power law, namely l/-\/t for T > and l/t when T = 0. Both behaviors are exponentially small in Te, and much 
smaller than subleading powers of Tg that were neglected already. So to the accuracy considered we can set M{t, H) 

equal to MeqiTe,H). 

When the field is slowly changed in the course of time, the leading behavior of /ii is 
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Ml 



(5.38) 



As long as |yUi| ^ HTe/K^ the main change of M(t) is still expressed by Meq{Te{t), H{t)). Since / ^ (1 + r)/t, this 
condition can be written as 



dH 



dint 



HKTejl + r) 
Y2 



dH 



dTe 



AHK{l + r) 



(5.39) 



This condition is re asona ble, and easily satisified near Tg — 0. 

When condition ( 5.38 ) is not fulfilled, it is not possible to describe U{t) and M{t) by Te{t) alone. One needs a 



second effective variable, namely the effective field He{t) = H{t) + SHe{t). Setting = y/f^ + H^, we then have 
from quasi-equilibrium formulae at {Tg, H^) 



u = —Kp 



2 



K 



H 
— ( 
K 



-5H„ 



2 



(5.40) 



m 



He 



(1 



2Ke 



H 

K 



5H, 



H 

2^2 



Tp 



From the definitions (5.12) we can make the identifications 



Ml 



^3 



6He 



or their inversion 



Te = 2e 



Mi; 



5He = - 



r2 -772 



Ml 



(5.41) 



(5.42) 



(5.43) 



One can now consider any class of fields that change logarithmically slowly in time. After solving the dynamics 
and Hp follow. One can then also study fluctuations and the fluctuation-dissipation relation for this more general 
case, and look for universal behaviors. This matter is the subject of currentresearch ||^, which falls outside the scope 
of the present paper. 



5. Fluctuations and fluctuation-dissipation relation 



again given by (4.52), 



The energy fluctuations still follow from (3.22), and depend on H only through K. They are, to leading order, 

(5.44) 



1 



TV 



This allows to solve the cross-fluctuations from ( |3.25| ), 



A 1 



1 + r-^ 



AK^ (l-r)(2-r) 



(5.45) 



Finally, up to corrections of order T^, the M-fluctuations satisfy an equation very similar to (4.44). This yields 



2K 



(5.46) 



where we notice that terms of relative order T^/A have cancelled on both sides of the equality sign. Considering the 
fluctuations in Mi = M, M2 = Mg, we have 



Cab(t, t) 



-H 



a-\-b—2 -p2T^ /J.^ 

A'3 ^ 2k' 



(5.47) 
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The difference with the oscillator model of previous section, is that all four terms (a, 6 = 1,2) now have corrections 
of relative order T^^/A^, that will be neglected from now on. They decay slower than the terms we keep, but we are 
not interested in the model-dependent very long-time regime, where they dominate. 

The time-dependence of the Cmm{t, t') follows from (3.34). In the interesting, not-very-asymptotic regime the term 
Cem can be neglected. This implies finally that 



Cab{t,t') = Cab{t', t') j^^J^ 

which involves h{t) defined as in eq. ( 4.39| ), with, very analogous to ( 3.19| ), 



fit) 



dxW{f3x) 



'2e{t){l-e/{2K) 
The equal-time correlators take the value 



p{x\e(t)) = d,Ap[Q\e){l + r){l 



2K 



2A 



Gab{t+,t) 



8A 



-H 



a+b~2 



piO\e) (1 



2A 



(5.48) 



(5.49) 



(5.50) 



Both Cab{t,t') and Gab{t,t') have a time dependence h{t)/h{t'). The fluctuation-dissipation relation ( p. 34 ) again 
holds with the same effective temperature (4.51) as in the oscillator model, 



Tejt) 
fit) 



2{T, - T) 
A 2T,-T 



(5.51) 



In deriving this result we noticed that terms of relative order Te/K, as appearing in eq. (5.47), have cancelled. We 
can now redo the consistency check of (2.35) and verify that, up to relative order T^, 



dt Gmmit,t ) — PeCmrnitjt) = ^7^(1 ^ 7777 ] 

K-^ IK 



Using that Te w A2iv:/(8 Int), we also find 



dH 



T.t 8K\nt if 2 



(5.52) 



(5.53) 



As mentioned before, this is one order of magnitude smaller than (5.52). The relation ( ^.33 ) now implies dm/dTf,\T,H 
~H/2K^, so that eq. ( 2.31) becomes 



fiuct ^ Q 

A mm ^ 



K^^ ^ 2K'^^^2I0'10 



(5.54) 



In view of the prediction ( 5.23 ), this is the desired answer to the considered order. 

Also the energy correlation- and response function are essentially the same as in the oscillator model. This i mplies 
in particular the fluctuation-dissipation relation for energy fluctuations of the non- universal form ( 4.55 ), ( 4.56 ). 

Notice that when there are no random fields, T = Q ^ K = H the energy and the magnetization are proportional 
to each other, viz. E = —HM. For comparison with realistic glassy systems, the model becomes too poor. The above 
formula for the magnetization seize to hold when F < Tg. For F — > one finds the magnetization correlations from 
the energy correlations. One then finds the relations (4.55), (f4.56|) both for energy and magnetization. 



6.. DISCUSSION 



In this paper we consider the question whether the glass transition can be phrased in a thermodynamic framework. 
In a series of letters we have given already several arguments in favor of this possibility j^] ]To[ |. The present, 
admittedly long, paper is meant to explain enough details of this approach to make the picture and its various steps 
and assumptions transparant. We do this by working out in detail two simple models, that, in our feeling, are closer 
to reality than mean field spin glasses. We have pointed out there that a minimal thermodynamic description needs 
one more parameter to describe the situation, which could be the age of the glass, or the cooling rate at which it was 
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formed. We have discussed that for thermodynamics a more useful variable is the effective or fictive temperature Te, 
introduced half a century ago by Tool . In this paper we notice that the basic result for the change of heat in a 
glassy system, 

dQ <TdScp + TedI, (6.1) 

immediately leads to a specific heat Cp — Ci + C2dTe/dT, that was Tool's starting point for the analysis in the glass 
formation region. In section IV. 1 we continue along his lines by studying in detail a certain non- linear cooling traject, 
proposed recently by us in ref. [|l^ . This cooling scheme is applicable to any glass forming substance and expected to 
give universal scaling curves of OTg/dT in the glass formation region, independent of the material considered, provided 
that the glass transition region is narrow. It would be most interesting to test this idea on a realistic glass forming 
liquid. One shou ld first determine, for once and for all, the equilibrium timescale Teq{T) and then do glass experiments 
of the type ( 4.14 ) within the considered range. It contains two parameters: the glass formation temperature Tg, where 



the cooling timescale becomes comparable to the equilibration timescale, and a "speed" parameter R. The resulting 
form for dT^/dT lies, after rescaling the width, on a universal scaling curve, that only depends on R. 

We have here worked out the situation where the effective temperature shows up as extra variable, though in 
principle it might be needed to consider as many effective parameters as there are macroscopic observables. We 
had already shortly considered the Ising chain with Glauber dynamics. In that model the non-equilibrium energy 
at zero-field can be described by introducing the effective temperature. Its definition then coincides by equating the 
Teq{Te) with t. The behavior at non-zero field appears to be non-universal 0. More or less the same happens in 
the backgammon model, for which the dynamics at zero field has been partly solved One could couple the 

system to a particle bath, and the chemical potential would play the role of an external field. So far this case remains 
to be worked out. We have, therefore, focussed on very simple, exactly solvabel models, namely the Bonilla-Padilla- 
Ritort model of Monte Carlo dynamics of uncoupled harmonic harmonic oscillators |^^, and our recent model of 
Monte Carlo dynamics of uncoupled spherical spins in a random field [ p^ . At low temperature both models have an 
Arrhenius law for the equilibrium timescale. Upon cooling from high temperatures, they will sooner or later fall out 
of equilibrium. Though it may come as a surprise for some readers, we have shown that these oversimplified models 
with their unrealistic dynamics still share in their off-equilibrium phase universal properties of realistic glasses and 
models for glasses. 

A description with only the effective temperature applies whenever the volume of the glass forming liquid (or the 
magnetization of the glassy magnet) is close to its quasi-equilibrium value set by the effective temperature obtained 
from fitting the energy. If this condition is fulfilled, the old objections against a thermodynamic description of the 
glassy phase can be inspected in detail. We have stressed that the most fundamental paradox, namely violation 
of the first Ehrenfest relation, is merely based on the misleading expectation that there is one, ideal value for the 
compressibility. This notion has arosen from equilibrium considerations for the glassy state. They donot apply by 
definition, and have hindered progress until our recent works in this field. Indeed, from the knowledge of spin glasses, or 
from the solution of the present models (where xzfc = 0), we know that the compressibility or the susceptibility can 
take a broad range of values immediately below the glassy transition. This means that no alternative determination is 
allowed, removing immediately the whole paradox: though half a century of research led to the general belief that the 
first Ehrenfest relation is always violated, it is actually satisfied automatically |^ . This point is underlined in figure 
2.1, where we present a 3d-plot of data for the glass transition in atactic polystyrene, collected in careful experiments 
by Rehage and Oels |^. Though these authors claim that the first Ehrenfest relation is violated, and then continue 
to investigate a modified version, we explain that it is satisfied. 

We have also pointed out that for glass forming liquids the Maxwell relation between dU /dp and dV/ dT is modified, 
which is not so surprising in view that equilibrium is not reached. The second Ehrenfest relation is then also modified, 
since it relies on the Maxwell relation. This fact implies that the Prigogine-Defay ratio can indeed be different from 
unity. We should recall that Davies and Jones showed that 11 > 1, while DiMarzio found that a deeper analysis 
of their equations leads to 11 = 1 Both approaches, however, are based on the assumption that at the glass 

transition an unspecified number of order parmeters freezes in, an assumption that was often made in the sixties and 
seventies. Such assumptions are invalid, however. What happens at the glass transition is that certain slow modes 
fall out of equilibrium, but on longer scales they may reach equilibrium again, even though other modes may then 
have fallen out of equilibrium. The upshot of this is that the Prigogine-Defay ratio can be different from unity. In 
contradiction to the standard belief, it can also be less than unity. We have pointed out before |^ that this already 
occurs in experiments on atactic polystyrene though this was long not recognized. 

It has been the important contribution of statistical physics to relate temporal fiuctuations in macroscopic observ- 
ables to their averages, the most known relation being C = < dH"^ >. It is natural to investigate whether such 
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relations have some universal-looking generalization in simple models for glasses. We have found that this is indeed 
the case for fluctuations in observables coupled to external fields see Section 2.5. These equations have been 
guessed with an eye on results from present models, in combination with some standard arguments on the short-time 
contributions. 

Our progress was initially hindered by the fact that such general formula appear not to hold for energy fluctuations. 
For the ferromagnetic Ising chain aging at T = from a random initial condition, we already realized that at zero field 
the energy defines an effective temperature Tg — 2Jln47rt, that coincides, to leading order, with the one following 
from equating the time scale with time, viz. Teq{Te) — t The energy correlations can be calculated from the 
Derrida-Zeitak spatial correlation function for the non-meeting of two random walkers on a line. Indeed, the Ising 
chain is mapped to random walking of interfaces by setting Si_i/2Si+i/2 = 1 — 2pi with pi = 1 if an interface is 
present at i, and zero else. It follows that U/N = J(— 1 -I- 2p) where p = 1/^ is the average density of interfaces, with 
^ — \/47rf = exp(2/3e J) the correlation length. It holds that 



l(<5L/2) = Aj'ip - + 8J2 ^(a,o - p') (6.2) 

j>0 

where Cy — (piPj) = C{rij) is the correlation function, given by eq. (59) with g = 2 in a paper by Derrida and Zeitak 



C{r) ^ p^{l-e-^^ +2ze-'' I e-""du), (6.3) 

where z = ^/tt pr. One thus finds 



-{SU')^Aj'p{3-2V2) (6.4) 

This differs from the naively expected quasi-equilibrium result T^dU /dTg ~ AJ^p by a numerical prefactor 3 — 2\/2 = 
0.17157288. We were informed by J.M. Luck that in the backgammon model a similar phenomenon occurs: the 
energy correlations are a factor two smaller than naively expected . In the models of the present work the energy 
fluctuations are even smaller by an order of magnitude in T^/A. These explicit examples show that there cannot exist 
a simple quasi- universal formula relating energy fluctuations with the specific heat. The underlying reason hereto is 
that in the energy the leading fiuctuations from different terms already cancel, leaving model-dependent, subleading 
effects only. 

After completion of the original manuscript, mr. L. Leuzzi verified that both in the harmonic oscillator model and 
in the spherical spin model there holds the following relation between the specific heat and the energy fluctuations: 



df 



1 ,..,9, dU 



AT, 



H,T V dT 



dU 



H,t J dT, 



dT, 



H,T \ dT 



H 



where T,^'^^ is given by eq. (4.56) for both models. The first two terms in the right hand side cancel, and so do, of 
course, the ±1 terms. In analogy with eq. (2.28), we could interpret this relation as C = (7^"'=* -|- (j^oss _|_ (^conf^ 



the factor 1/T in the first term could imply that this attempted generalization is special to our present oscillator and 
spin models, as is already expressed by the fact that T,^'^'^ T,. It was also realized that for the cross derivatives 
dU/dH and dM/dT similar relations either do not occur, or are much more complicated. The same holds for the 
fluctuation-dissipation relations connected to these two quantities. 

Nevertheless, fluctuations in observables that couple to global external flelds, appear to behave in a universal 
way, at least to leading and dominant subleading order. These fluctuations are interrelated with the off-equilibrium 
fluctuation-dissipation relation (FDR). Originally observed by Horner [ pl] - fi3[ an d then extended by Cugliandolo and 
Kurchan , this has become a popular test for glassiness of model systems [|42|] [ pT| p2| . For the models consid- 
ered the effective temperature showing up in the FDR is essentially the same as the one occurring in thermodynamics. 

Absence of a universal relation between energy fluctuations and the specific heat is very welcome for gravitating 
systems, which often have a negative specific heat. We indeed showed that the present approach can immediately be 
applied to phrase the laws of black hole dynamics in a non-equilibrium thermodynamic framework p8[ | . The role of 
bath temperature is played by the cosmic backgound temperature, that of the effective temperature by the Hawking 
temperature, that of the configurational entropy by Bekenstein's black hole entropy, while the short times processes 
have no sizeable entropy. These ingredients bring a very close analogy with the picture discussed here, and are not 
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even disturbed by the fact that the specific heat has the "wrong" sign. Let us mention, however, that negative specific 
heats are in no way hmited to gravitation: they also occur in the present, extensive sohd state models, when one 
heats up the system in the glassy phase, as expressed by the negative values of dT^/dT = 2C in figures 4.1-4.4. 

Our picture for thermodynamics of the glassy state thus connects macroscopic observables via the first and second 
law, and relates their derivatives with respect to external fields with their fluctuations, thereby embedding the FDR- 
effective temperature in a larger thermodynamic framework. It is expected to be valid for a yet unknown class of 
glassy systems. Let us close this discussion by mentioning that very recently numerical data in the glassy phase of a 
binary Lennard- Jones system were interpreted in terms of an effective temperature, that dominates the (short-time) 
vibrational properties, in full harmony with the picture proposed above H]. 
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